Geometrical interpretation of the Casimir effect 
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Casimir forces are a manifestation of the change in the zero-point energy of the vacuum caused 
by the insertion of boundaries. We show how the Casimir force can be efficiently computed by 
consideration of the vacuum fluctuations that are suppressed by the boundaries, and rederive the 
scalar Casimir effects for a series of the Dirichlet geometries. For the planar case a finite universal 
force is automatically found. Consistent with other calculations of the effect, for curved geometries 
divergent (non-universal) expressions are encountered. They are interpreted geometrically following 
Candelas and Deutsch (1979) as largely due to the divergent self-energy of the boundary contributing 
to the force. This viewpoint is supported by explicit calculations for a wedge-circular arc geometry in 
two dimensions where non-universal and universal contributions into the effect can be unambiguously 
separated. We also give a heuristic derivation of the purely geometrical expression (Sen, 1981) for 
the non-universal piece of the Casimir energy due to an arbitrary smooth two-dimensional Dirichlet 
boundary of a compact region. 

PACS numbers: 03.70.+k, ll.10.-z, ll.10.Gh, 42.50.Pq 



I. INTRODUCTION 

Casimir interactions are due to the macroscopic re- 
sponse of the physical vacuum to the introduction of 
boundaries. They were first derived as an attractive 
force between perfectly conductive parallel plates induced 
by the zero-point motion of the electromagnetic field [l[ 
There is convincing experimental evidence for the reality 
of these forces Q and a vast body of literature dedicated 
to various aspects of the phenomenon [3| . 

The Casimir interaction £ is the difference between the 
vacuum energy of the system constrained by the bound- 
aries and that of free space. Since boundaries made of 
real materials are transparent to sufficiently high-energy 
modes [lj, the high energy spectrum is unaffected by 
the geometry of the system, and only a finite range of 
the spectrum need be considered. However, in the theo- 
retical treatments of this effect the vacuum energies are 
usually calculated from an effective low-energy harmonic 
field theory (such as quantum electrodynamics in the case 
of the electromagnetic Casimir effect), so that they are 
approximated by the sum of zero-point energies of a col- 
lection of simple harmonic oscillators with a spectrum 
u) = c|k| (where c is the speed of light). In this model, 
the dispersion relation holds for arbitrarily large wave 
vectors k; both the "constrained" and "free" vacuum en- 
ergy densities are ultraviolet divergent; and the Casimir 
interaction is the difference between two infinite quanti- 
ties. The theory resolves this problem by a soft-cutoff 
modification of the large-k part of the spectrum that 
leads to a finite vacuum energy. The result is not very 
sensitive to the form of the cutoff, so that the Casimir 
interaction can then be extracted by taking the cutoff 
to infinity at the end of calculation. Other approaches 
to the calculation that make use of analytic continuation 
[^ and dimensional regularization [5| techniques give the 
same answer, thus adding to the credibility of the result. 

As already noted, the divergences are more mathemat- 



ical artifacts than physical reality. The important virtue 
of the model is that for many geometries a finite univer- 
sal result depending only on h, c, and macroscopic length 
scales is obtained without introducing a cutoff. However, 
this is not always the case: specifically, the divergences 
occurring for spherical geometry in even space dimen- 
sions do not cancel [a 0]. What this means physically 
represents an open problem; it seems to imply that in 
two dimensions a conducting ring placed in vacuum is 
unstable. 

The goal of this paper is two-fold: first, we show how 
the Casimir effect can be efficiently computed by direct 
consideration of the fluctuation modes that are elimi- 
nated by the presence of the boundaries. Technically 
this is accomplished by using a method closely related to 
the path- integral approach to the Casimir effect [8| • 

Second, we compute the Casimir effect in a two- 
dimensional wedge-circular arc geometry. The results 
of our analysis lend support to the idea of Candelas 
and Deutsch [9| that the divergences encountered in the 
case of curved boundaries are of geometrical nature and 
mostly due to divergent self-energy of the boundary con- 
tributing to the Casimir force. Since the geometry we 
study is a relative of the ring geometry in two dimen- 
sions, our results also shed some light on the physical 
meaning of the divergent Casimir force exerted on the 
ring in two dimensions [g, [7| • 

The organization of this paper is as follows. In Section 
I we introduce a method of computation of the Casimir 
effect which is based on consideration of the vacuum fluc- 
tuations eliminated by introduction of boundaries. The 
efficiency of this technique is illustrated in Section II 
where the scalar Casimir energies and pressures are re- 
derived for a series of standard geometries: planar ge- 
ometry in an arbitrary dimension, circular geometry in 
two dimensions, cylindrical and spherical geometries in 
three dimensions, and finally, spherical geometry in an 
arbitrary space dimension. We observe that the Casimir 



pressures are always divergent in curved geometries. Fol- 
lowing Candelas and Deutsch Q we pursue an interpre- 
tation of these divergences as non-universal contributions 
due to the geometry of the boundary. This is made ex- 
plicit in detailed calculations of Section III where the 
effect is computed for the wedge-circular arc geometry 
in two spatial dimensions with the wedge edges subject 
to the Dirichlet and periodic boundary conditions. The 
latter includes the circular geometry in two dimensions 
as a special case. Our main observation here is that non- 
universal contributions to the Casimir effect can indeed 
be attributed to the geometry of the boundary. Moreover 
the non-universal pieces of the effect turn out to be insen- 
sitive to the change of the topology of the vacuum while 
the universal ones are. This fact is used in Section IV to 
heuristically re-derive a purely geometrical formula [10( 
for the non-universal piece of the Casimir energy of an ar- 
bitrary smooth Dirichlet boundary of a compact region. 
We conclude (Section V) by summarizing our findings. 



II. FORMALISM 

In what follows we will be analyzing Casimir effects in 
a Gaussian field theory with the Euclidian action 



S[w] 



h l T f dm 

dTd d x( C - 2 (^) 2 +(Vw) 2 



(1) 



where T is the temperature. The real scalar w is a func- 
tion of the d-dimensional position vector r and imagi- 
nary time r, and is periodic on the Matsubara circle, 
u>(r,0) = w(r,H/T) pTJ . The action ^ is applicable at 
energies low compared to some scale Hluo; the frequency 
cutoff loq is provided by the properties of the material 
the boundary is made of. 

The scalar field theory (TTJ) can be viewed as a toy ver- 
sion of electricity and magnetism. The divergences en- 
countered for curved geometries exist in both theories, so 
they are not due to specifics of electricity and magnetism. 
Therefore we study the problem in the simpler setting of 
the scalar theory. 

The zero-point energy can be calculated by means of 
a functional integral that makes use of the correspon- 
dence between the Feynman path integral for the d- 
dimensional field theory and the partition function for 
a d + I-dimensional classical statistical mechanics prob- 
lem. The path integral is 



Z,, 



Dw(r 7 T) exp(—S[w]/h) 



(2) 



but it can be interpreted as the partition function for a 
classical statistical mechanics problem with the Hamil- 
tonian S at a fictitious temperature which is equal to 
Planck's constant [12j. The zero-point energy corre- 
sponds to the "free energy" per unit "length" in the imag- 
inary time direction, so that £q = —h(lriZ w )/(h/T) = 
— TlnZ,,,. 



Now assume that the vacuum is disturbed by the ad- 
dition of sharp boundaries which constrain the field in 
some way and thus eliminate some degrees of freedom of 
the vacuum fluctuations. The constrained field (which we 
will refer to as v) inherits any boundary conditions im- 
posed on w as well as new conditions on surfaces Di and 
Nj of Dirichlet (v\i) i = 0) or Neumann (dv/dn\N. = 0) 
type, respectively, where the subscripts i and j label the 
boundaries and dv/dn is the normal derivative. We will 
write the difference between the original and constrained 
fields in the form w(r, t) = v(r, t) + u(r, t), where u(r, t) 
can be chosen to satisfy the d + I-dimensional equation 



(■ 



d 2 



' 2 dr 2 



-A)u = 0, u\ Di =fi(r,T),—\ Nj =g j (r,r) (3) 



where fi and <?j are functions defined on the boundaries 
and determined by the boundary values of w; they play 
the role of dynamical variables of our approach. The 
reason for defining u this way is that it eliminates the 
cross term in the action, so that S[w] = S[v] + S[u]. Then 
the functional integral factors into integrations over v and 
u, so that Z w — Z V Z U . Since the zero-point energy for the 
unconstrained system is determined by Z w and the zero- 
point energy for the constrained system is determined 
by Z v , the Casimir energy is given by £ — T\nZ u . This 
result allows us to calculate directly the change in energy 
due to the fluctuation modes that have been eliminated. 
A simplification is achieved by expanding all the dy- 
namical variables of the problem into a Fourier series 
in the imaginary time domain; for example u(r, r) = 
J^ w u^ (r) exp icuT where the Fourier coefficients u u (r) 
are solutions to the boundary-value problem for the 
Helmholtz equation 



(A--^-)w^ = 0, u u \ Di =/o,,i(r), -q^\n, =gu,j(r) (4) 



The calculation of the action S[u] is further simplified 
when the identity (Vn) 2 = div(uX7u) — uAu is substi- 
tuted into Eq.((T|). Then the integral of div(uVu) over 
d d x transforms into a sum of surface integrals. The re- 
maining integral over dr vanishes due to the relation 
Au = —d 2 u/c 2 dr 2 and the condition of periodicity, 
w(r, 0) = w(r, h/T). As a result we find 

1 r h / T r 

S[u] = - / rf^2_/ { u ^ u ]idsi 

= 2T 51 / K Vli -^]* ds » ( 5 ) 

uj.i 

Here [4>] i stands for the discontinuity of ip across the i-th 
boundary, and the summation is performed over all the 
boundaries. 

Although our approach is applicable to an arbitrary 
number of the Dirichlet and/or Neumann boundaries, in 
all the cases considered in this paper only one Dirichlet 
boundary contributes into the integral (J5]). Since we are 



dealing with a harmonic field theory, the solution to the 
boundary-value problem ^ is linear in the surface field 
fa, and thus the action Se is a quadratic diagonal form 
of f u - 



Se -^^2 



\u\ 



2T^A„(M/c) 



(6) 



where the subscript v (representing one or more indices) 
labels the normal modes of the field u that have been 
eliminated by the boundary in question, and Aj,(|w|/c) > 
is a set of geometry-dependent coefficients. Then ac- 
cording to the prescription outlined above, the Casimir 
energy is given by 

p T^ x ^TK{\u\/c) 



2 ^ 

He yr-^ 

2^^ 

v 

he y^ 

~2^^ 



dftln \ v (k) 



kcLkF{ — )— — In A,, 
luq an 



(7) 



where in the second step we took the zero-temperature 
limit according to the rule ^~J W — > (h/T) J du>/2ir and 
dropped contributions that are independent of geometry. 
In the third representation, a monotonic cutoff function 
F(y) rapidly decaying for y > 1 and satisfying F(Q) = 1 
is made explicit; the cutoff prescription, hen/2 = Hui/2 — > 
(hui / 2) F (uj / u>o) , reflects penetrability of the boundary to 
high-energy modes [10J] . In most of the formulas below the 
cutoff function is suppressed and invoked only as needed; 
analysis of such cases is conducted for an arbitrary cutoff 
function. We note that the last representation of hints 
at a relationship of our approach to a contour integral 
method of calculation of the Casimir energy [3] ■ 



A. Planar geometry 

Consider three Dirichlet planes at z — 0, z = a, and 
z = L, where z is one of the axes of the d-dimensional 
rectangular coordinate system and < a < L. We are 
interested in the Casimir pressure exerted on the mid- 
dle partition at z — a. The outer boundaries are fixed 
in place, so that there is no need to look beyond them. 
Since the space is uniform relative to translations paral- 
lel to the boundaries, the field w w (r) is expanded into a 
Fourier series u u (r) = ^ q w w q(z) expiqr^ where rj_ is 
the position vector perpendicular to the z axis. Then the 
boundary-value problem ((4| for the Fourier coefficients 
u^q(z) becomes 

d 2 uj 2 
(-7-2 -Q 2 2")«uiq = °' u ^qlo,i = °> Uwqla = fuq (8) 



The particular solution to (jSJ) is 



Uuq(z) = f t 



^q 



sinh(|«;|z) 
sinh(|ft;|a) ' 



^ z ^ a, k 2 — q 2 H 5- 



U) 

~ 2 



, . „ sinh(|«;|(L — z)) r ,„. 

«u,q(«) = fvq . , \) T ((, a < z < L 9 

smh(|«|(L — a)) 

Substituting this in Eq. ([5]) we see that only the partition 
at z = a contributes to the action Se with the result 

Se = ^H |«|(coth(|«|a) + coth(|«|(£ - a)))|/. q | 2 
w,q 

(10) 
where A is the macroscopic (d— l)-dimensional area of the 
boundary. Eq. (|10[) has the required form Eq.([B]), with co- 
efficients that become large (A becomes small) for large n, 
because for high frequencies or large transverse wavevec- 
tor the disturbance introduced by the fluctuating bound- 
ary condition is localized at the boundary, to within a 
length that is proportional to A itself. This will lead in 
what follows to a divergent surface energy. 

The Casimir energy per unit area can be deduced from 
Eq. (fl7J|) according to the rule Eq.([7j), with the result 



£ T v^, /coth(|/da) + coth(k|(L-a)) 

7 = ^L ln 



A 2A 



2-kT/HA\k\ 



uj.q 

~* ~\ J d \2n) dq ln ( coth (N«) + cothd^Ki. - a))) 

= — / K d ~ 1 dKhx(coth(Ka) +coth(K(L - a))) 

2 Jo 

(11) 

where in taking the macroscopic limit we used the rule 
Eq ~> Ajd d ^ 1 q/(2Tr) d ^ 1 . Additionally, in going from 
the first to the second representation we took the T = 
limit and dropped all a-independent contributions not 
influencing the pressure on the partition. Hereafter all 
such contributions will be systematically dropped; how- 
ever, we should note that in the present case part of 
what has been omitted is a divergent integral represent- 
ing an infinite surface energy. The parameter Kd in the 
third representation is the surface area of a d-dimensional 
unit sphere, S d = 2n d / 2 /T(d/2), divided by (27r) d . The 
Casimir pressure on the boundary, V = —d(£/A)/da can 
be found in closed form 



V = 



(47t)t 



he 



1 



(L-a) 



d+l 



-,d+l 



(12) 



where T(x) and C( x ) are Euler's and Riemann's gamma 
and zeta functions, respectively [131 ]. In arriving at 
(I12p we used the gamma function duplication formula, 
T(z)T(z + i) = 2 1 - 2z ^T(2z) [H and the value of the 
integral [14f 



a; d (cothx - l)dx = 2" d r(d + l)((d + 1) (13) 



We see that the partition at z = a is attracted 
to the closest outer boundary; specifically, in one di- 
mension Eq. (fT2")) reduces to the well-known result [ljj. 



Since the outer boundaries impose the Dirichlet bound- 
ary conditions we can imagine joining them together. 
Then Eq. (fTS)) describes Casimir interaction between two 
boundaries; taking the L — > oo limit we then reproduce 
the result of Ambj0rn and Wolfram [5j. 



B. Circular geometry 

Consider a Dirichlet circle of radius a in two spatial 
dimensions. The boundary- value problem (|U) for this 
geometry becomes 



ld_ _d_ 1 8 2 
p dp dp p 2 dtp 2 



= 0,Uu;\a = fu(<P) 



(14) 

where p and <p are the polar coordinates. Seeking the par- 
ticular solution in the form, u u (p,(p) — i? w (p) expin<p, 
where n is an arbitrary integer, we find that the radial 
function R^(p) satisfies the equation 



d 2 R ul 1 dRu 



dp 2 p dp 



(=V + \)R* = o 



(15) 



whose linearly-independent solutions are modified Bessel 
functions I n (\ui\p/c) and K n {\ui\p/c) [13]. Thus the par- 
ticular solution to the boundary- value problem (TTJ} finite 
at p = and decaying as p — ► oo is 



C. Cylindrical geometry 

The Casimir interaction for three-dimensional in- 
finitely long Dirichlet cylinder of radius a can be in- 
ferred from corresponding results for the circle. First, 
expand the field u u (r) into a Fourier series, u u (p, <p, z) = 
J2q u ^qz(Pi ( P) e:>c P^ ( lz z where p, tp, and z are cylindrical 
coordinates, and the z axis coincides with that of the 
cylinder. The Fourier coefficients u^ qz (p, ip) are solu- 
tions to the boundary value problem (JT4J) except that the 
role of ijj 2 jc? is played by q 2 + uj 2 /c 2 . Thus cylindrical 
analogs of Eqs.(fr5|)-(rr7]) can be written down by replac- 
ing the summation index lu by the combination (to,q z ), 



and a; /c by k 2 



2 /c 2 . Moreover, an extra fac- 



tor of the macroscopic cylinder length £ appears in the 
counterpart to Eq. fTTj) . As a result the Casimir energy 
per unit length £/£ and the pressure on the boundary 
V = — (2Tra£)~ 1 d£/da are given by 



£ 

£ 



, oo 

he y^ 

4^ ^ 

72= — OO 



/"OO 

/ kcIkIii (I n {K,a)K n {hLa)) (20) 

/o 



, OO />OC J 

V = ~^h £ / ^dx^{\n{I n {x)K n {x))) (21) 

n— — oo ^ 

respectively Again, these results are known [17| ■ 



y^ In(\u\p/c) 

^ K n (\w\p/c 
£- K n {\u\a/c 



fun exp irvp, p^a 



■f un expin<p, p> a (16) 



n— — oo 



At the circle p = a both of these reduce to fuif) — 
E^L-oo /w« ex P i nL P- Substituting the solution (fTr?)) in 
Eq.([S]) and performing the angular integration we find 



airh sr-^ 



i f 4(M a / c ) -MM a /c) \ | , ,2 

l \I n (\oj\a/c) K n {\Lo\a/c)j Uun] 



(17) 
We observe that for large n or large ui the prefactor in this 
sum is large, because the functions Eq. ((T6|) are localized 
near the circle boundary. This result can be simplified by 
application of the Wronskian relationship K n (z)I n (z) — 
I n (z)K n (z) = 1/z 13]. The Casimir energy £ and the 
pressure V = — (2iTa)^ 1 d£/da are then implied by Eq.j7|) 



. OO 

2?r ^ 



n— — oo 



dnln (I n (Ka)K n (Ka)) 



(18) 



f- °° />0O 7 

i^3 E 7 xdx-MIn(x)K n ( X )) (19) 



n— — oo 



The results (fT5]) and (TIT))) are due to Sen [T(| (see also 
Refs. [|,[3). 



T8. 2 8. 1 
(r — ) + — 



D. Spherical geometry in three dimensions 

Consider a Dirichlet sphere of radius a in three spatial 
dimensions. The boundary- value problem Q for this 
geometry becomes 

\_d_ . d_ _L 8_ 

r? V sin 9 89 {Sm 88* + sin 2 6 dip* 

uj 2 -\ 
- -Z juu, = 0, u,„| = /„(«, yj) (22) 

where r, 9 and </? are spherical coordinates. Seek- 
ing the particular solution in the form, u u (r, 6, <p) — 
R,.>(f)Yi m (9,ip) where Yi m (6,tp) is a spherical function 
of order I, m [la] , we find that the radial function R u (r) 
satisfies the equation 



d R^ 
dr 2 



2 dR UJ 
r dr 



1(1 + 1) 



Ru =0 



(23) 



are 
Thus 



whose linearly-independent solutions 

/, + i(Mr/c)/VF and tf I+ i(Mr/c)/Vr 0. 
the particular solution to the boundary-value problem 
finite at the origin r = and decaying as r — > oo is 



oo I 



VEl l+i (\u\r/c) 
Z^ Z^ ~tz~t 1\. A „i^ JulrnYl m {9,(p), r sC a 



1=0 m=-l 

oo l 



y/fl l+ i{\u\a/cy 



2^ Z^ ^"^^ — n. A _i_-, JuimYi m {v,<p), r>a 



1=0 m =-l 



y/rK l+ i(\uj\a/cy 



(24) 



At the sphere surface r — a both of these reduce to lem ((4]) that has to be solved is 



fw(0,V>) = X)j )TO /a>Zm*im(0>¥>)- Substituting the solu- 
tion l[2"4")) in Eq.([5]) and performing angular integration 
we find 



S F = 



ha 



2T ^ 771 



\f« 



i(\u\a/c)K l+ i(\u\a/c) 



(25) 



In arriving at (|25[) . similar to the circular case, we used 
the Wronskian of the I t+1 / 2 (z) and K t+1 / 2 (z) pair as well 
as the property of orthogonality of the spherical functions 
[l§ |. As a result the Casimir energy £ and the pressure 
on the boundary V — — (Aiva 2 )~ 1 d£ /da are implied by 
Eq.dU) 



z=o 



1) 



dfiTn 



/, + i(Ka)if,,i(Ka) 



(26) 



V 



Pi _ _ Z* 00 ^7 

^J2(2l + 1) xdx-ln 



d , Ij+ 1(^)^+1(3;) 



(27) 

These equations reproduce the results due to Bender and 
Milton |6| who also derived more general relationships for 
a Dirichlet sphere in space of arbitrary dimension. 



E. Spherical geometry in arbitrary number of 
dimensions 



As a last demonstration of the efficiency of our tech- 
nique we derive the expressions for the Casimir energy 
and pressure for a Dirichlet sphere of radius a in ar- 
bitrary space dimension. Before proceeding we remind 
the reader [2fJ, [2l|, l22j that in d dimensions spheri- 
cal functions can be arrived at by first enumerating all 
homogeneous linearly-independent harmonic polynomi- 
als of degree I, P/^( r )- After the transformation to 
spherical coordinates X\ = rcost^i, x 2 = rsin^x cosip 2 y 
..., x d -i = rsin^isin^ 2 ---sin( < 9 d _2COS( < 9 d _i, x d = 
rsanpi sin^--- sin ^-2 sm(p<i-i, the polynomials take 
the form Pi,d{r) = r l Yi t d(9), where Y^d is a (/-dimensional 
spherical harmonic of degree I and 9 refers to the point 
on unit sphere with angular coordinates <pi, (p 2 ,...,ipd-i- 
For / ^ there exist 



mi,d 



(2l + d-2)(d + l-3)l 
(d-2)!Z! 



(28) 



linearly- independent spherical harmonics of order I. Fol- 

r(m). 



lowing Mikhlin [22j, these will be denoted by Y^ (9), 

m = l,2,...,mi t d- For an y I the functions Y^'(ff) can 
be made orthonormal on the unit sphere, and we assume 
this is the case. 

The remaining steps mirror our discussion of the two- 
and three-dimensional cases. The boundary-value prob- 



92 d-ld 1 . uj 2 \ 

TTJ H TT + -^&6 ? )U U = 0, U u \a = Jlu(V) 

or 1 r or r z c z J 

where Ag stands for a multiplicative angular piece of 
the Laplacian; explicit expressions for A# in two and 
three dimensions are given in Eqs. (fT4]) and (|22|) . re- 
spectively. Seeking the particular solution in the form, 
u u (r, 9) = R w (r)Y l ^ (9), we find that the radial function 
Ru}{r) satisfies the equation 



d 2 R u} d-ldR^ 



dr 2 



dr 



lu 2 1(1 + d-2) 



R^ = (30) 



where we employed the result AgY"; ™ (9) = —1(1 + d — 
tyYi™ (®) [20, [2l|. Linearly-independent solutions to 
([3H)l are r~^~ I l+ d-2_(\uj\r/c) and r^^~K l , d-2 (\u)\r/c) 
[23j thus implying that the particular solution to the 
boundary-value problem (I29|) finite at the origin r = 
and decaying as r — ► 00 is 



EE 

2=0 m=l 

00 mi,d 

EE 



i =0m =i r 2 Ii + i^(Ma/c) 
00 m w a^rK 



i+- 



,(|«|r/c) 



■( m )^( m ) 






(|w|o/c) 



(31) 



At the sphere surface these reduce to fu>(9) — 

Sz m fujiuT (^) which is a series expansion in spherical 
harmonics of an arbitrary function defined on a sphere 
[2TI . HI]. Substituting the solution J3T]) in Eq.© and 
performing the integration we find 



St 



ha d - 2 ^ 

2T V Ijjl*-' 

1,) I rn * 1 ■-) 



17 



( m )|2 
I I 



(\u)\a/c)K l+ d=2(\u>\a/c) 



(32) 



As a result the Casimir energy £ and the pressure 
on the boundary V = —(Sda d ~ 1 )~ 1 d£/da (here Sd — 
2it d l 2 /T(d/2) is the surface area of a unit sphere) are 
given by 



£ 



hr J2. r°° Ii 1 d-2 (na)K, , d-2 (na) 

T^'-'L dKl " (1^ (33) 



v 



he 



2nSda 



d+l 2-^1 



mid 



=n 



r°° d . Ii+^i(x)K l+ ^2(x) 

x / xdx— In ?^ ^5 — 2 (34) 

/o 



dx 



a; 



which are the results due to Bender and Milton [G] ; specif- 
ically the expression for the pressure is identical to their 
Eq.(3.5). 



In three dimensions Eqs. ([33[) and ([34]) clearly reduce 
to Eqs. PS]) and (|27p. In order to see that the two- 
dimensional results are also reproduced, we recall that 
22] for d — 2 and I > there are only two linearly- 
independent harmonic polynomials of degree I, namely 
the real and imaginary parts of {x\ + 1X2) 1 : while for 
d = 2 and I = there is only one polynomial which 
is a constant. This implies that in two dimensions the 
degeneracy factor ([28]) satisfies the rules, to,;. 2 = 2 for 
I > and too, 2 = 1- With this in mind we see that for 
d = 2 Eqs. fll§ and $M$ reproduce Eqs.CEB]) and CEH]). 



III. DIVERGENCES AND THEIR 
GEOMETRICAL INTERPRETATION 

Although Eq. (JT2J) gives a finite universal Casimir pres- 
sure for the planar version of the problem, its circular 
(fl"9]) . cylindrical {JSTJ) , spherical ([27]), ([34]) (and any curvi- 
linear) counterparts are divergent. Various techniques 
have been successfully used to remove the divergences in 
the case of a cylinder [17| and in spherical [6( geome- 
tries of odd dimensionality predicting universal Casimir 
pressures. However no consistent removal procedure was 
found capable of handling the spherical geometry of even 
space dimension [6] . This includes the experimentally rel- 
evant two-dimensional case where as an alternative, Sen 
[101 ] proposed to view the action |T]) as an effective low- 
energy theory to be supplemented by a cutoff function 
like in Eq.©. This removes the divergence, and leads to 
a finite non- universal effect. 

In order to understand the difference between the pla- 
nar and curved geometries we notice that the Casimir 
force is the change of the energy upon infinitesimal 
displacement of the boundary. A geometrically sharp 
boundary possesses a divergent energy per unit area, 
coming from the exclusion of the high frequency and 
short wavelength modes from an increasingly narrow re- 
gion near the surface. This divergent self-energy does 
not contribute to the Casimir force in the planar case 
because the overall area remains fixed as the boundary is 
displaced. However this is not the case for curved bound- 
aries. Indeed a change of the radius of a circle implies 
a change of the perimeter and as a result the divergent 
self-energy will contribute into the force. This idea orig- 
inally due to Deutsch and CandelasQ was recently re- 
expressed by Graham, Jaffe and co-workers |2J] and by 
Barton [25| . The implication is that a real- world curved 
boundary may be responsible for large non-universal por- 
tion of the physically measurable Casimir force [26j . This 
is consistent with Sen's observation [Hj] that the coef- 
ficients of the non-universal terms contributing into the 
Casimir energy are geometrical objects such as the length 
of the circular boundary, etc. We also note that in the 
high-temperature limit the Casimir free energy can be 
shown to have entirely geometrical nature - it is expressed 
in terms of a surface integral of a quadratic function of 
local curvature 1271. 




a 



Figure 1: Wedge of opening angle /3 with superimposed arc 
of radius a in two dimensions. 



The geometry-dependent contributions to the Casimir 
effect can be classified into three groups: (i) non- 
universal terms, which diverge as a power of the cutoff 
frequency u>q, (ii) nearly- universal terms, which diverge 
as a logarithm of loq; and (iii) universal terms, which 
remain finite as luq — > 00. The nearly-universal contri- 
butions are characterized by universal amplitudes and a 
weak dependence on the choice of the cutoff function. 

A further test of the idea of Refs.(9l.[24l.[25t would con- 
sist in calculation of the effect in a curved geometry where 
various contributions into the Casimir interaction can be 
unambiguously separated. A natural candidate is a sys- 
tem characterized by more than one macroscopic length 
scale where the area or circumference of the boundary is 
independent of its curvature. 



A. Wedge-circular arc geometry in two dimensions; 
Dirichlet case 

Consider a wedge of opening angle (3 with superim- 
posed arc of radius a, Fig. 1. This is a geometry in 
which both a Casimir force and torque exist. 

The semi-circular {§_= 7r) version of this problem has 
been studied earlier [28| using the zeta-function regular- 
ization technique. Similar to the circular case, a diver- 
gence was found which led the authors to conclude that 
" for obtaining the physical result an additional renormal- 
ization is needed" . 

The wedge-arc configuration with an arbitrary angle (3 
has been considered in Ref . 29] with the aim of revealing 
the regularities in the boundary non-smoothness contri- 
butions to the heat kernel coefficients. The local charac- 
teristics of the vacuum have been discussed by Sakharian 
and collaborators for a scalar field with Dirichlet bound- 
ary condition in general space-time dimension [301 ] and 
for the electromagnetic field in three-dimensional space 

a 

Although the Casimir interaction can be obtained from 
the Casimir energy density of Ref. [301 ] , our approach pro- 



vides a quicker route to the final result. Indeed the 
Casimir energy can be inferred from the results for the 
circular geometry since the boundary-value problem we 
need to solve is closely related to (fT¥|) . The difference 
is that we seek a solution inside the Dirichlet wedge 
< ip < (3, thus implying u UJ {p 1 Lp) = R UJ (p)sm(Trrup/ (3), 
n = 1,2, ... for the particular solution. The radial func- 
tion Ruj(p) satisfies the same Eq. flTS)) with n being re- 
placed by im/ (3. As a result the solution in question can 
be obtained from Eq. (flU|) by replacing the order n of the 
Bessel functions with nn/ ' f3, the angular function exp imp 
with sin(Tmip / P) , and restricting the summation over n 
from unity to infinity. The calculation of the Casimir en- 
ergy is similar to that for the circular geometry with the 
result 
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where in the second representation we restored the cutoff 
function (see Eq.((7])) and introduced the length of the arc 
s = (3a. We then employ the uniform asymptotic n»l 
expansion of Debyc 13,32]: 

2n(l + yY /2 I„(ny)K n (ny) = l + ^{ I l 

6 5 -| _L_ 

(1 + y 2 ) 2 + (1 + y 2 ) 3 / + [ n 3 > 

which can be used to evaluate the energy 
j3 <C 1 limit. To leading order we find 
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where in going from the second to the third representa- 
tion we used Eulcr's sum 
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Although the expression (|37|) is divergent in the wo -^ oo 
limit, it is important to realize that it depends on a and 
(3 only through the arc length s — {3a. In order to inter- 
pret (|37p geometrically we note that the dimensionless 
combination ujqs/c entering the argument of the cutoff 
function is a ratio of the macroscopic s and microscopic 
c/loq length scales of the problem, and naturally one has 
s 3> c/ujq which is the range of applicability of the the- 
ory. Then in order to evaluate the first term of (|37p we 



expand the cutoff function in a Taylor series around zero. 
The leading F(0) = 1 term of the expansion then gives 
an integral of the form (|13| , thus generating a universal 
7r/ic/48s contribution to the Casimir energy. All higher 
order contributions of the expansion are convergent and 
proportional to negative powers of the cutoff frequency - 
in the ojqs/c 3> 1 limit they all vanish. The second term 
of Eq. (f3"T)) can be ignored because it leads to a finite 
geometry-independent energy contribution of the order 
hw 0l while the contribution from the last term is linear 
in the arc length. As a result we arrive at 
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where a is the line tension coefficient; as expected, it is of 
the order of the microscopic energy scale hwo divided by 
the microscopic length scale c/wo ■ If we specialize to the 
case of an exponential cutoff function (F(t) = exp(— t)), 
then our expression for the line tension reproduces Sen's 
result a = ^ 2 /(4ttc) [13]. 

Even though Eq. (f3"9"| is the leading order term of the 
/3« 1 expansion of the Casimir energy, there is no re- 
striction on the segment length s = f3a since the radius a 
can be arbitrary. We see that to the leading order there 
is only the dependence on the segment length (thus im- 
plying locality) but not on its curvature. Eq. ((39|) also 
predicts that if the system were given the freedom to 
choose the optimal value of s by minimizing the energy, 
this would be a microscopic s ~ c/uq which is a com- 
promise between the energy decrease due to shrinking of 
the segment length (the non-universal term) and repul- 
sive universal piece of the Casimir interaction preferring 
large segment lengths. The precise value of the optimal 
s implied by Q39p is not to be trusted as s ~ c/luq is at 
the verge of applicability of the theory. If we take the 
simultaneous (3 — > 0, a — > oo limit but keep the length 
s = j3a fixed, the geometry of Fig.l turns into that of 
an infinite strip of width s whose sides are connected by 
a straight Dirichlet bridge; its Casimir energy is given 
by Eq. (|3"9")) . From this viewpoint the universal 7rfi.c/48s 
piece of (|39p can be understood by approximating the 
bridge by a strictly one-dimensional Dirichlet interval. 
The Casimir energy of a free field confined to an interval 
of length s is — 7rftc/24s and attractive [la ]. In our case 
it is repulsive as these fluctuations are eliminated by the 
boundary. Additionally, the magnitude of our effect is 
smaller since the one-dimensional approximation incor- 
rectly assumes an infinitely sharp field localization and 
thus a larger change of the vacuum energy. 

In order to improve on Eq. (|39|) we subtract £v-> , 
Eq.((37D, from Eq.©. The outcome U = £ - £ (1) is 
expected to be universal or nearly-universal since the 
largest non-universal fraction of the effect is already in- 
cluded in £ W: 
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where we suppressed the cutoff function. Using the ex- 
pansion (|36[) again we evaluate the energy (T4T))) to the 
next order in (3 <C 1. If the integration is performed first, 
then we find 



U = - 



256tts ^-J 
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(41) 



which is marginally divergent. This means that setting 
the cutoff function at unity is not justified, and that there 
is a weak dependence on the cutoff frequency wq. There- 
fore we restore the cutoff function, go beyond the leading 
order term in the Debye expansion l]36[) and write Eq. ([40]) 
as 
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The sums over n can be computed with the help of Euler's 
formula (|38[) and two consequent relationships: 
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where the prime stands for the derivative with respect to 
x. Comparing the right- and left-hand sides of Eqs.(j38|). 

([43]) . and flS} we observe that the sum J2 {••■} in Eq. ([42"l) 
approaches a constant limit as x — > 0. This implies 
that the integral over x is convergent at the lower limit. 
On the other hand, in the x — > 00 limit we find that 

Y^ \ ■■■ \ —* — 1/(2 7 £) with the implication that if the cut- 
off function is set to unity, then, consistent with Eq. (j4"T1l . 
the integral over x would logarithmically diverge at the 
upper limit. The role of the cutoff function consists in 
restoring the convergence by effectively setting the up- 
per limit of the integration at ujqs/c 3> 1. In view of 
the logarithmic character of the divergence, the outcome 
is rather insensitive to the upper integration limit. As 
a result we arrive at the nearly-universal contribution to 
the Casimir energy 

25D7TS c 

which has logarithmic accuracy. 

The result (|45|) can be also understood heuristically 
starting from the divergent expression (|4ip and employ- 
ing the equality Yl n =x n ~ X ~ C + In iV + e/v where C is 
Euler's constant and e^ — ► as N — > 00 [13| . The param- 
eter N is estimated by recalling that the suppressed cutoff 
function F(irnct/uJos) effectively ends the sum of 1/n at 



n = N such as irnct/cuos « 1 and the relevant value of 
t, as implied by the second-order term of the Debye ex- 
pansion (|36j) . is of the order unity. Thus N ss luqs/c S> 1 
and we recover Eq. (|4"5|) . 

The result (|45|) implies that for fixed arc length s when 
the leading order contribution (|55|) can be regarded as 
a constant, the zero-point motion induces a widening 
torque —dU/d[3 which could be expected on physical 
grounds. Combining Eqs. ([39|) and (|45|) we arrive at the 
expression for the Casimir energy of the Dirichlet wedge 
of arc length s and opening angle (3 = s/a <^ 1 

hes cj s irhc 

t = as ;r in (46 

2567ra 2 c 48s K ' 

This result supports a geometrical interpretation for 
the non-universal contributions to the Casimir energy. 
Indeed, when the frequency cutoff u>o tends to infin- 
ity, the strongest (oj 2 .) divergence comes from the first 
term of (|4"6")) proportional to the length of the bound- 
ary thus supporting the idea of Refs.[3, [2J, |25| , The 
sub-leading logarithmically divergent term of (|46[) can 
be also understood geometrically as being proportional 
both to the length of the boundary and to the square 
of its curvature. This term can be alternatively inter- 
preted as a finite-size correction to the surface tension: 
a — > a — hcln(ujQS / c) /256na 2 . We note that Eq. pS)) has 
a local character and exhibits a term-by-term correspon- 
dence with its circular counterpart [101 ] - 

The accuracy of the Debye expansion (f3"6"| . implies that 
Eq. ptj]) approximately captures the whole ^ f3 ^ 2ir 
range. 



B. Wedge-circular arc geometry in two dimensions: 
periodic boundary conditions 

The Casimir effect is sensitive to changes in space 
topology. In our next example we demonstrate that the 
non-universal contributions to the Casimir effect do not 
seem to be sensitive to topology while the universal ones 
are. Consider the wedge-circular arc geometry as in Fig. 
1 but assume that the wedge edges are physically iden- 
tical. This is implemented by imposing the condition of 
periodicity, u{p, ip) — u{p, if + /3), where the opening an- 
gle (3 can now be arbitrary. This means that the field u 
belongs to a conical surface with a Dirichlet circle which 
is a distance a away from cone's apex as measured along 
the cone surface [2{|. For [3 = 2tt the geometry becomes 
that of a plane with the Dirichlet circle. 

The Casimir energy of this configuration can be in- 
ferred from the solutions of the circular and wedge cases. 
Indeed, the boundary-value problem we need to solve 
is posed by Eq. (fT4"l) with the boundary field satisfying 
the condition of periodicity, f^i^p) = f^i^P + /?)• This 
implies u ul (p 1 (p) = Ruj(p) sin(27m</5 / '(3) , n = 1,2,... or 
u^lp, ip) = Ruj{p) cos(2irrup / (3) , n — 0,1,2,... for the 
particular solution. In both cases the radial function 
Rui{p) satisfies the same Eq. (Tr5t with n being replaced 



by 2irn/ (3. The Casimir energy can then be inferred from 
Eq.® as 
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where the first term is due to the fluctuating angle- 
independent Fourier component of the boundary field 
fu while the sum is a contribution from the angle- 
dependent components. The latter are represented by a 
sine and a cosine Fourier series contributing equally into 
the Casimir energy. We note that for (3 = 2tt Eq. ([4"T)) 
reduces to its circular counterpart, Eq. (fT8|) . 

In order to understand the geometry dependence of the 
first term of (j47|) we restore the cutoff function F(y): 
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The derivative removes the dependence on the opening 
angle /?. Since 2xIo(x)K (x) ss 1 [13| for large values of 
the argument, the integral <|48 j) diverges linearly at the 
upper limit when the cutoff function is set at unity. The 
presence of the cutoff function effectively sets the upper 
integration limit at aj$a/ 'c thus leading to a geometry- 
independent estimate £q ~ fuvo- The a-dependent part 
can be found by subtracting from (|48|) the same integral 
with Iq(x)Ko(x) replaced with its large argument limit 
l/2x: 
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This is the universal piece of the Casimir effect when the 
cutoff function is set at unity (i. e. the luocl/c — > oo limit 
is taken). The integral was numerically evaluated by Sen 
[irj | and more accurately by Milton and Ng [16[ whose 
value is used in (|4"9"|) . 

The second term of Eq. ([47)) can be analyzed in a 
manner similar to that of Eg. ([55)) . The results of such 
analysis combined with (|49p give the Casmir energy for 
the wedge-arc geometry with periodic boundary condi- 
tions: 
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where the coefficient of surface tension is defined in 
Eq. ([3"9")) . We see that the cutoff-dependent terms of 
Eas. (|46| and ([50]) are the same. At the same time the 
effect of topological change in going from the Dirich- 
lct to periodic boundary conditions manifests itself in 
the universal parts of the Casimir effect. Indeed, the 
s-dependent piece of the interaction, irhc/12s, is four 
times larger than its Dirichlet counterpart, 7rfic/48s, in 
Eqs. ([3"9")) and (|4"o"|) . This mirrors the relationship between 
the Dirichlet and periodic Casimir energies for the free 



field in one dimension [33|. Additionally, the last attrac- 
tive term of ([50]) is unique to the periodic geometry. 

For the special case of a circle in a plane, s = 2wa, Eq. 
([5TI| simplifies to 



£ = 2-Kua - 
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This agrees with Sen's result [10( except for the magni- 
tude of the last term. We believe Sen's value 0.045 is a 
typographical error (according to the numerical results 
he quotes, the value should be 0.054). 

Similar to the argument following Eq. (|39|) we can take 
the simultaneous [3 — > 0, a — > oo limit but keep the 
segment length s = (3a fixed. Then Eq. ([50"]) simplifies to 
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which is a periodic counterpart of Eq. (|39p : it gives the 
Casimir energy of a Gaussian field defined on an infinite 
cylinder of radius s/27r in the presence of a Dirichlet circle 
belonging to the cylinder surface. 



IV. ARBITRARY SMOOTH BOUNDARY OF A 
COMPACT REGION 

Since for the circular geometry the length of the bound- 
ary and its curvature are not independent, a route differ- 
ent from direct calculation is needed to demonstrate the 
geometrical nature of the non-universal contributions to 
the Casimir interaction, Eq. (|5ip . For that purpose Sen 
[lOJ I employed a result due to Pleijel [34j which relates 
asymptotic behavior of a certain integral of a Green's 
function to geometry. In the process Sen derived a re- 
markable expression for the non-universal part of the 
Casimir energy £ (u>o) for an arbitrary smooth boundary 
of a compact region. Here we show how this result can 
be understood heuristically by employing our findings for 
the wedge geometry. 

First we notice that for the wedge whose edges are 
Dirichlet/periodically connected by an innnitesimaily 
small arc/circle of length ds the non-unversal part of the 
Casimir energy (see Eqs. (f4*T|) . pp| and (|5"0"|l ) can be 
rewritten as a differential relationship 



d£(u>o) — ads + dU = ads 
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where C — 1/a is the curvature of the arc. Our calcula- 
tions indicate that the energy ([53)) has to be attributed 
to the arc itself rather than to the conditions imposed 
on the wedge edges; those conditions determine the uni- 
versal part of the Casimir effect. With this in mind let 
us consider an arbitrary smooth boundary of a compact 
region. Then Eq. (|53p can be interpreted as representing 
a non-universal contribution into the Casimir energy due 
to the infinitesimally small boundary element ds whose 
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position, the length s along the boundary, is measured 
relative to an arbitrary reference point and C(s) is the 
curvature at location s. 

For an arbitrary smooth boundary the non-universal 
piece of the Casimir energy can be found by integrating 
([53]) over the length of the boundary 



calculations indicate that additivity does not hold for the 
universal parts of the Casimir energy which could not be 
expressed in a differential form. The implication is that 
it would be very difficult if not impossible to come up 
with a purely geometrical formula for the universal part 
of the Casimir effect. 
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where S — J ds is the total length of the boundary. We 
notice that the marginally divergent sum of n _1 does not 
have any reference to specific geometry. This leads us to a 
conjecture that (similar to the case of the circular arc) the 
effect of the frequency cutoff ojq consists in replacement 
of the divergent sum in ([54]) with the finite logarithm of 
the large ratio of the macroscopic and microscopic length 
scales: 
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Since Eq. (|55[) has logarithmic accuracy, the precise mean- 
ing of the macroscopic length scale is not essential - we 
have chosen it to be the total length S of the boundary. 
The result (|55p is due to Sen [jjj, [35| . Our heuristic argu- 
ment which leads to Sen's result ([55)) also links the Pleijel 
formula [3J| to the Debye expansion ([36]) which is in the 
heart of our wedge results. 

The remarkable feature of Eci. ([55j) is that it is solely 
determined by local geometry of the boundary. In ar- 
riving at (|55[) we employed additivity of the interaction 
which is known generally not to be the case for Casimir 
interactions. We note however that the additivity only 
holds for the strictly non- universal piece of the effect. For 
the nearly- universal second term of Eq. ([55]) the additiv- 
ity is logarithmically weakly violated. This allowed us 
to proceed by assuming the additivity and to deal with 
the logarithmic divergence separately. Additionally, our 



V. SUMMARY 

To summarize, we have demonstrated how Casimir ef- 
fects caused by sharp boundaries can be efficiently com- 
puted by focusing on the quantum fluctuations elimi- 
nated by these boundaries. The applicability of this 
method is not limited to the scalar field theory (fT]), 
Dirichlet boundaries, zero-temperature limit or to the ge- 
ometries we have considered. 

Second, we presented an explicit calculation to support 
the idea [3, [24], [25| that the divergent Casimir forces en- 
countered in the presence of curved boundaries have ge- 
ometrical origin - they are largely due to divergent self- 
energy contributions. Our analysis also supports Sen's 
proposal 10] that the Casimir effect in a Dirichlet ring in 
two dimensions is finite and non-universal with the cutoff 
frequency o>o supplied by the properties of the material 
the boundary is made of. More work is needed to fur- 
ther explore the geometrical nature of the non-universal 
parts of the Casimir effect as they may be responsible 
for largest contributions into experimentally measurable 
Casimir force. 
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